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We report application of a method for extracting gravitational waves to three-dimensional 
numerical simulation on coalescing binary neutron stars. We found the extracted wave form 
includes the componets corresponding to the quadrupole part in the Newtonian potential of 
the background metric, if it is monitored at a position not far from the central stars. We 
present how to eliminate it. 



We are constructing computer codes on three-dimensional numerical relativ- 
> ■ ity.^)'2) At first we used the conformal slicing condition, in which the metric be- 

oo : comes the Schwarzschild one in the outer vacuum region so that if the three-metric 

^ I is split into the Schwarzschild background and the perturbed parts, the latter can 

■ be considered as the gravitational waves at the wave zone.'^^ However, it was found 

\ that this slicing involves unstable modes and long-term evolution of coalescing bi- 

' nary neutron star cannot be followed. ^^'^^ Then we started to construct a new code 

using the maximal slicing condition. In this slicing, the perturbed part of the three- 
metric includes gauge dependent modes and therefore we need gauge-invariant wave 
extraction. Recently gauge-invariant wave extraction methods have been given as 
nonspherical perturbations of Schwarzschild geometry. In this letter, we report 
0^ \ application of a method based on them to three-dimensional general relativistic sim- 

' ulation on coalescing binary neutron stars. 

We use (3-|-l)-formalism of the Einstein equation and write the line element as 
^ \ ds"^ = —a^dt^ + ^ij{dx^ + (3^dt){dx^ + l3^dt). Outside of the star, we split the total 

' spacetime metric into a Schwarzschild background and non-spherical perturba- 

■ - - ' tion parts: g^^, = g\^^ + h^^} + h^^}, where g^i^^ is a spherically symmetric metric 

given by 

gl^Uxf^dx" = -N'^dt^ + A^dr"^ + R^ide"^ + sin^ ^#2) (1) 
and h^^} and h^^^} are even-parity and odd-parity metric perturbations, respectively; 
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43 = "^^ Sin^ 9 [KimYim + Glm {YlmM " Wlm)] , (3) 
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-h^til^Yim,^/ sin 9 
-h"il^Yim,4^/ sine 
h^L^im/sine 
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\ sym sym 



h^"^ Y, 



YimfisinO \ 
sinO 



sym 



-h^L^ira sine 



(4) 



(e) (e) 

where the symbol 'sym' indicates the symmetric components, Hn^n, h^i^, hi/^, Ki^, 



Gim, h^Qiin^ ^^'^ functions of t and r; Y/^ is the spherical harmonics, 

Xi^ and Wim are given by 

Xlm = '^{Ylm,e4> - Ylm,<t> COt 9) , Wir, 



Yim,de - Yim,d cot 9 - Yim,,j>4,/ sin^ 9. (5) 



Prom the linearized theory about perturbations of the Schwarzschild spacetime, the 
gauge invariant quantities and •P'^^^ are given by^^ 



and 



'^itV, r) = -^2{A-2)/A {ArN^k2ira + Arkum)/{A + 1 - ^N^) 



for the odd and even parity modes, respectively, where A = l(l + 1), 



(6) 
(7) 

(8) 

(9) 



The quantities !P'*-°'* and ^Z'*-^-' satisfy the Regge- Wheeler and the Zerilli equations, 
respectively.^'^) Two independent polarizations of gravitational waves ^+ and hx are 
given by 

^+-i^x =-^B'^2(^>0 + '^tH*>0)-2iln., (10) 

l,m 



kiim = Kim + N\Gim,r - 2N^h^li>Jr 



and 



where 



lYim = {Wim - iXirn/ sin e)/^/A{A 



(11) 



In numerical calculations, the functions N'^{t,r), ^^(t, r) and B?{t,r) of the back- 
ground metric are calculated by performing the following integration over a two- 
sphere of radius r:^) 
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-)dr2, (12) 



where df2 = sin 9d9d(l). The components of metric perturbations are 

H2lmit,r) = ^ J grrYCmdn, 



(13) 
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where * denotes the complex conjugate. 

We need angular integrals over spheres for constant r, such as 



Im 



df?, 



(14) 
(15) 
(16) 
(17) 

(18) 



r=ro 



f{x,y,z)df2= / /(ro,e,(/))df?. 



(19) 



If numerical simulation is performed using Cartesian coordinate system, we need 
interpolation to obtain the values of f{rQ,9,<j)) from f{x,y,z) at the grid points. It 
is, however, not easy to fully parallelize the procedure on a parallel computer with 
distributed memory. We therefore rewrite Ea. H19() as a volume integral, namely. 



-^(^o) = / f{x, y, z)6{r — ro) d^x = lim 



1 



V^rar^ 



F{x,y,z)e-'^'-'''>^'/^' d^x, (20) 



where r = y^r^ + y^ + z^. Numerical integral with a = Ax/ 2 gives a good value to 
Ea. (|2U() . where Ax is the separation between grid points. 

We performed numerical simulation for a coalescing binary consisting of two 
identical neutron stars of mass 1.5Mq and evaluated the gravitational waves. The 
details of our code will be shown elsewhere^^-* but it is essentially the same as Refs. 2) 
and 12). The lapse function and the shift vector are determined by the maximal 
slicing and the pseudo-minimal distortion conditions, respectively. We used uniform 
475 X 475 X 238 Cartesian grid with Ax = IMq assuming the symmetry with respect 
to the equatorial plane. As for an equation of state, we use the 7 = 2 polytropic 
equation of state. The initial rotational velocity is given so that the circulation of 
the system vanishes. The ADM mass of the system is 2.8Mq. 

Figure^shows the gravitational wave forms r/i_|_ and rhx on the z-axis evaluated 
at r = 110, 120, 130 and 14OM0 as functions of the retarded time t — r. The lines 
of rhjf- x{t — r) estimated at r = 110 ~ 14OM0 for t — r > coincide with each 
other. Then the waves proportional to and propagating at the speed of light are 
extracted. For i — t < 0, however, h includes a non-wave mode proportional to 
as shown in Fig. ^ where r^h is plotted as a function of t. This mode corresponds 
to the quadrupole part in the Newtonian potential of the background metric. It 
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Fig. 1. Plots rh+^x along z-axis at r = 110 ~ 
14OM0 as a function of t — r. 
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Fig. 2. Plots r'^h+^x as a function of t. 



decreases fast as the merger of stars proceeds. Since this mode is proportional to 
while the wave mode is to r~^, the former will be negligible if the waves are 
monitored at a few times farther position. Now we can eliminate the non-wave mode 
from and using Fourier transformation as follows: 

• Assuming that total waves are expressed as a sum of wave parts F(t — r)/r and 
non-wave parts G{t)/r^, 

• Fourier components of h{t, r) are written as 

Kir) = -^F^{r) + ^G^{r). (22) 

where 

F^ = ^ J F{t)e-'^'dt and = ^ J G{t)e-'''' dt. (23) 

• From the values of h^{r) in different radial coordinates ri and r2, can be 
given by 

p ^ 4Kir2) - rlKjri) 

• By inverse Fourier transformation, we can get the gravitational waves, which 
do not include non-wave modes, 

^—F^e"^' duj. (25) 

The resultant wave form is shown in FigEl The curves represent the average 
of /14. and hx calculated at r = 110, 120, • • • 2OOM0 and twice the dispersion 2a is 
shown as error bars. 
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Fig. 3. Wave forms rh^^x along z-axis as a func- 
tion of t~r. The curves are averages of rh+^x 
estimated at r — 110 ~ 2OOi\f0 and error 
bars denote 2a. 
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Fig. 4. The comparison /i+,x defined by 
Ea. ll25ll and ft+,x defined by Ea. ll26ll ob- 
tained with the conformal sHcing (c) and 
the maximal slicing (m). 



Here we define h+ and as 
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and /i> 



(26) 



respectively, wliere h 
and d) 



= (P lij - Oij 

(det(7jj))T2 . The pseudo- 
minimal distortion condition demand- 
ing dt{djhij) = guarantees /i+^x to be 
transverse-traceless if djhij = at t = 0. 
It is our case since we assumed the ini- 
tial three-metric to be conformal flat, 
hij = 0. Then they can be considered 
as the gravitational waves on z-axis in 
the conformal slicing, while they include 
gauge dependent modes in the maximal 
slicing.^-' To compare /i+,x with /i+,x 
for the conformal slicing as well as the 
maximal slicing, we perform numerical 
simulation for a coalescing binary of two 
M = I.OM0 neutron stars. As shown in 
Fig. 12 they almost coincide with each 
other, while a small deviation is found 
in /i+,x in the maximal slicing. Then 
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Fig. 5. The energy spectrum of the gravita- 
tional waves plotted in Fig El The curves 
are averages of AE/Aui estimated at r = 
110 ~ 2OOM0 and error bars denote 2a. 



we found that the gauge mode in /i+^x 
is small even in the maximal slicing. 

Finally to investigate a possibility that the excitation of the quasi-normal modes 
can be seen by the numerically calculated waves, we evaluated the energy spectrum 
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(27) 



where is the Fourier transformation of ^^^{t^r). Figure [S] shows the energy 

spectrum of the waves plotted in Fig. |31 The fundamental frequency of / = 2 for 
a Schwarzschild black hole of mass 2.SMq is u; = 25 msec~^. A peak near this 
frequency appears in Fig. |SJ Unfortunately, however, rotating angular frequency 
just when the merger of the stars finishes is 12 ~ 15 msec~^ and thus they will 
radiate the waves of frequency near lj = 25 msec~^. So that more precise calculation 
is necessary to discuss whether this peak corresponds to the emission of the quasi- 
normal mode of the formed black hole. 
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